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COMPUTATION  OP  THE  RADIATION  PATTERN  OF  A  REFLECTOR  ANTENNA 
INTRODUCTION 


The  radiation  pattern  of  a  reflector  antenna  is  a  function  of  the  distri¬ 
bution  of  the  surface  current  on  the  reflector.  Once  this  surface  current  is 
known,  the  antenna  pattern  can  be  computed  directly  by  integrating  this  cur¬ 
rent  distribution  over  the  reflector  surface.  This  surface  current  is  induced 
by  the  incident  field  from  the  primary  feed  source.  If  the  primary  feed 
pattern  is  known,  it  is  straight  forward  to  compute  the  incident  field  on  the 
reflector  surface  and  hence  the  induced  surface  current.  This  process  in 
general  involves  an  integration  of  the  surface  current.  In  the  past,  the  pri¬ 
mary  feed  is  located  at  the  focal  point  and  its  radiation  pattern  is  known 
in  terms  of  the  reflector  coordinate.  The  induced  surface  current  on  the  re¬ 
flector  is  hence  readily  to  be  written  down.  However,  in  some  cases,  the  pri¬ 
mary  feed  is  not  necessarily  located  at  the  focal  point  and  its  pattern  may  be 
known  at  a  different  coordinate  system.  It  is  thus  necessary  to  have  a  coor¬ 
dinate  transformation  before  the  reflector  surface  current  can  be  determined. 

Another  aspect  of  this  problem  is  that  the  Green’s  function  of  this  inte¬ 
gration  involves  both  the  coordinates  of  the  surface  current  and  the  coordi¬ 
nates  of  the  field  points.  Each  time  the  radiation  field  of  a  different  field 
point  is  computed,  this  integration  must  be  performed  over  again.  Unfortu¬ 
nately  this  integration  in  general  cannot  be  integrated  into  a  closed  form. 
Repeating  of  this  numerical  integration  is  time  consuming,  particularly  when 
a  large  number  of  field  points  must  be  computed.  In  this  report,  a  method  is 
presented  in  which  the  Green’s  function  is  expanded  into  an  orthogonal  series 
which  can  be  Integrated  into  a  closed  form.  The  secondary  radiation  pattern 
is  the  summation  of  this  series.  Hence,  it  requires  only  to  perform  this 
integration  once  to  determine  the  coefficients  of  the  series.  It  thus  greatly 
reduces  the  required  computation  time. 

RADIATION  FIELD  OF  A  REFLECTOR  ANTENNA 

The  radiation  pattern  of  a  reflector  antenna  is  determined  by  the  surface 
current  distribution  on  the  reflector.  The  radiation  E  field  is  related  to 
the  surface  current  by  the  following  relations: 


.-Jko  *//& 


exp  Ok0P*R)da 


e”Jko* ffr  *  exp  Uk0P*R>d. 


Where  R  is  Che  distance  from  Che  focal  point  (reference  point)  of  the 
reflector  to  the  far  field  point  and  ?  is  a  vector  pointing  from^the  refer¬ 
ence  point  to  a  point  on  the  reflector  (Figure  1) .  While  o  and  4  are  re- 
spectively  the  unit  vector  of  the  spherical  coordinate  of  the  reflector, 
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equations  (la)  and  (lb)  can  be  written 


%  -  ®  • t 


(2a) 


i 


$  -  i 

-  j  S 


(2b) 

(2c) 


It  is  thus  straightforward  to  compute  the  radiation  pattern  of  a  reflec¬ 
tor  once  the  surface  current  distribution  is  known*  This  surface  current  is 
induced  by  the  incident  field  which  can  be  found  from  the  following  relation 


?s  -  2  (nxH± )  -  2^|  [nxCp'xEj)]  (3) 

o 

*  + 
where  n  is  the  unit  vector  normal  to  the  reflector  surface  and^p '  is  a  radial 

vector  from  the  primary  source  to  the  reflector  surface  while  E.  is  the  inci¬ 
dent  E-  field. 


Let  us  assume  that  the  primary  feed  is  far  away  from  the  reflector  sur¬ 
face  and  the  field  incident  on  the  reflected  surface  can  be  considered  at  a 
far-zone  region.  Furthermore,  mutual  coupling  effect  can  be  neglected.  Under 
these  conditions,  the  incident  field  on  the  reflector  can  be  represented  in 
the  following  form.  [1,2] 


where  P  is  the  total  radiating  power  from  the  primary  feed  and  G.(6',  <&')  Is 
the  radiation  pattern  of  the  primary  feed  which  usually  is  known  in  terms  of 
the  feed  coordinator  9'  and  The  quantity  p'  is  the  distance  from  primary 

feed  to  the  reflector  surface.  Vector  e.  is  the  unit  vector  which  defines  the 
direction  of  this  Incident  feed  on  th£  reflected  surface.  If  the  primary  feed 
is  linearly  polarized  at  a  direction  y,  this  vector  e.  is  then  in  the  direc¬ 
tion  u  and  normal  to  vector  p'.  This  can  be  expressed 


e±  -  p'x(yxp') 


Inserting  equation  (4)  into  equation  (3) ,  one  finds  that 


/cf(e;*’) 


exp  (Jk^p  * ) 


where  C 


(5) 


(6) 
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where  0q  is  the  limit  angle  of  the  reflector  dish. 

If  the  reflector  is  a  paraboloid,  its  surface  can  be  described  as 

X2  +  y2  -  4f (f-z)  (8) 

The  origin  of  this  rectangular  coordinate  is  at  the  focal  point,  and 


This  n  is  defined  in  the  unprimed  or  reflector  coordinates.  However,  vector 
and  ei  are  in  the  primed  coordinates.  Furthermore,  the  integral  is  per¬ 
formed  in  terms  of  the  0  and  <j>  coordinates,  while  the  G-  function  is  given 
in  terms  of  0f  and  <J>f.  Therefore  a  means  of  transforming  from  the  primed 
coordinate  to  the  unprlmed  coordinate  and  a  means  to  convert  0V  to  0  and  <J>' 
to  <p  must  be  provided.  In  the  next  section,  we  shall  discuss  the  coordinate 
transformation. 

COORDINATE  TRANSFORMATION 

The  coordinate  transformation  can  be  performed  most  conveniently  in 
rectangular  coordinates.  Let  us  assume  that  the  origin  of  thg  primed  coordi¬ 
nates  is  displaced  from  the  unprimed  coordinates  by  a  vector  £,  then  in  terms 
of  rectangular  coordinates  any  point  specified  in  the  primed  coordinate 
system  can  be  represented  in  the  unprimed  coordinate  system  by  the  relation 


3 


« 


where  £  ,  and  e  and  £  are  the  components  of  vector  £  defined  in  the  un¬ 
primed  coordinate*  z 


Besides  this  translation,  the  coordinate  system  may  also  be  rotated* 

This  rotation  can  be  conveniently  characterized  by  three  Euler  [3,4,5]  angles. 
The  transformation  has  the  following  relation. 


One  notices  that  e  ,  e  ,  £  are  defined  in  terms  of  unprimed  coordinates* 
By  use  of  the  above  refat itn»  Iny  point  defined  in  either  coordinate  system 
can  be  transferred  to  other  coordinates ,  provided  the  Euler  angles  and  e  are 
known. 
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When  the  Integral  I  of  Eq.  (7)  is  calculated,  the  integration  is  per¬ 
formed  in  the  unprimed  coordinate,  however,  the  primary  feed  pattern 
function  G^(0',  4>f)  is  defined  in  the  primed  coordinates.  Furthermore,  the 
distance  from  primary  feed  to  the  reflector  surface  current  point  to  be  inte¬ 
grated  (p')  is  also  in  the  primed  coordinates.  Therefore,  for  a  given  point 
on  the  reflector  surface  with  known  0,  <J>  and  p,  one  must  find  the  correspond¬ 
ing  0 1 ,  <j>’  and  p'  in  the  primed  coordinates.  This  can  be  converted  as 
follows.  For  a  given  point  with  known,  p,  4>  and  <J>,  it  can  be  represented  in 
the  rectangular  coordinate  as 

p  *  psin0cos<{>x  +  psin0sin<j>y  +  pcosGz  (14) 

This  can  be  transformed  into  primed  coordinates  by  use  of  Eq.  (13),  the  p  is 
then 

P  =  fjS'  +  f2yf  +  f32f  (15) 


where 

fi 

psin0cos<}>-ex 

f2 

-  [A]'1 

psinOsin^-e^ 

f3 

pcos  -e 

z 

and  xf,  yf  and  zf  are  unit  vectors  in  the  primed  coordinate.  Meanwhile  this 
same  point  can  be  represented  in  parameters  0f,  <{>f,  and  p*  as 

p  *  p fsin0 fcos<J>fxf  +  p'sinG'sin^'y* +  p'cosG’z1 


Hence  one  has 

p’sinQ’cosiJ)1 
pfsin0f sin$f 
p  fcos0f 


fx(0,4),p) 

f2(e,*,p) 

f3(0.^,p) 


(16a) 

(16b) 

(16c) 


One  can  solve  for  0f,  <t>’  and  pf  by  use  of  these  three  simultaneous  equations. 
The  second  problem  involved  in  this  integral  is  the  cross  product  of  the 
vectors  n,  and  e.«  Both  0*  and  e.  are  defined  in  the  prime  coordinate, 
but  n  is  defined  in  the  unprimed  coordinate.  Hence  the  cross  product  of  (S' 
and  e'  must  be  first  represented  in  a  rectangular  coordinate  of  the  primed 
system  which  are  then  transferred  into  the  unprimed  coordinate  system  by  use 
of  Eq.  (12).  After  this  transformation  is  completed  one  finds 


r 


-jWli  f  o  r 

C  exp(-jkoR)Jo  Jo 


[f  (e;*)x+f  (e,*)y+f  (e.cul 

x  y 


expfjk  p'  (e,(J>)  J  .  -  2  e 

- 2 - -  •  exptjk  p  *R]p  Sin6s«c  (17) 


5 


Since  the  vector  I  is  defined  in  the  rectangular  coordinate  system  which  has 
the  same  coordinates  of  the  field  point  P,  it  is  straightforward  to  compute 
the  components  Eg  and  E^. 

Our  task  is  then  to  evaluate  an  integrand  of  the  form 


exp[jk  p 1 (8,$) ] 

f(0,*)  - 

p' (6,*) 


(18) 


^  2  0 
exp{jkoP*R]p  sinecos  j  <39<i4> 

Effect  Of  Off-set  Feed: 

If  the  feed  location  deviates  from  the  focal  point  is  small,  pf  can  be 
approximated, 

p*~  p  -  p •  t 

^  p  _  e  sin0cos<j>  -  e  sins  in*  -  e  cosfl  (19) 

—  x  y  z 

•+  A 

Insert  this  into  equation  (18),  and  expand  the  dot  product  p  •  R,  one 
finds 


iSltH  exp(jL)  rd0d$ 

P» 

where  r  «  p  sin  0 

L  -  -k  (pcosdcos 0  +  p  +  e  cose)  +  k r  sin  0  cospcos* 
o  z  o 

+k  rsin0  sin$sin$  +  k  e  sin9cos$>  +  k  c  sinBsintf 
o  ox  o  y 

For  a  paraboloid, 


2f 

p  *  l-cos8 


Equation  (21)  can  be  manipulated  into  the  following  form. 


(20) 

(21) 

(21a) 

(21b) 

(21c) 
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~2fk  *  k  pcos0(l-cos6)  +  k  rsin9cos<J>cos$  +  k  rsin0sin<f>sin$ 
o  o  o  o 

+korevC0S<*1  +  korevSin*1  +  k_re'cos0cos4i  +  k  re  'cos0sin(Ji 

w  «  w  y  ox  o  y 

+k  e'  2fcos0 
o  z 


Ex  =  Ex/2f 

(23a) 

ey  =  Ey/2f 

(23b) 

where 

t '  *  e 

2  z/2f 

(23c) 

Let 

sin  6  'cos*'  *  sin  6  cos$  +  c' 

X 

(24a) 

sin  0  'sin**  =  sin  0  sin$  +  c' 

y 

(24b) 

These 

two  equations  represent  the  effect  of  off-focal  point  field  which 

shifts 

the  radiation  pattern  from  a  point  0  and  $  to  a  point 0f  and  . 

Sec 

L  *  pcos0 (cos0  -  1) 

(25a) 

[  2  2 

M  *  V  e  +  e 
v  x  y 

(25b) 

n  -  tan-1  e*  /€*x 

2 

(25c) 

^  2f*^  1  2f 

L  =  r  - - — 5-  (Mcos(4>-n)-c  —  ) 

c  *  z  r 

i+(^> 

(25d) 

Insert 

these  relations  into  equation  (2) f  one  finds 

r  ^  2it  ,* 

I  *  k|  °  |  F(r,<{>)exp[  jk  L  ]exp[jk  rsinO'cos 1 )  ] 

Jo  Jo  o  a  o 

Where 

rdrd$ 

(26) 

F(r,*)  *  (  exp[-j2k  f]expt-jk  Lj 

0  O  C 

(27) 

P' 
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If  the  primary  feed  is  located  at  the  focal  point,  then  e  ,  e  ,  and  e  are 
all  zero.  Under  this  condition,  0’  -  0,  -  <t>  and  L  *0.x  Clearly  thi  effect 

of  off-focal  point  feed  not  only  introduces  additions?  phase  delay,  it  also 
changes  the  beam  point  direction. 

Integration  of  equation  (26)  involves  both  0,$  and  0,  $  variables.  This 
appears  in  two  places.  The  first  place  is  the  exp[+jk  L  ]  term.  For  this 
term  one  may  calculate  this  integration  at  a  given  0  angfe  for  all  points  at 
the  vicinity  of  this  angle.  Larger  errors  will  certainly  appear  at  points 
which  are  remote  from  this  angle.  However,  there  is  a  way  to  correct  this 
error.  This  will  be  discussed  later.  At  the  time  being,  we  shall  assume 
that  is  computer  at  a  given  0  and  it  is  only  a  function  of  0. 


The  second  place  is  the  exp[jk  rsinO’cos^-^')  ]  term  which  we  shall  show 
in  the  next  section  that  it  can  be  integrated  into  a  closed  form,  if  F(r,<J>) 
is  expanded  into  an  appropriate  orthogonal  series. 

One  notices  that  the  integral  of  equation  (26)  involves  variable  0  and 
Thus  for  each  field  point  computed,  this  integration  must  be  repeated  all 
over  again.  For  simulations  of  a  reflector  antenna  or  computing  the  radiat¬ 
ing  pattern,  it  usually  involves  a  large  number  of  field  points.  The  re¬ 
quired  computer  time  is  astronomical.  To  overcome  this  difficulty,  Galindo- 
Israel  and  Mittra  [6]  proposed  to  expand  F(0,<{>)  into  an  orthogonal  series  so 
that  equation  (27)  can  be  integrated  into  closed  forms  as  a  function  of  0 
and  <$  of  the  field  coordinates.  Since  the  series  is  orthogonal,  the  coeffi¬ 
cient  of  the  expanded  series  can  be  easily  evaluated.  Galindo-Israel  and 
Mittra  showed  that  this  series  converges  very  fast,  only  a  few  terms  are 
needed.  This  will  greatly  reduce  the  required  computer  time.  This  proce¬ 
dure  will  be  summarized  as  follows: 

Let's  set 

r  -  sr 

o 

where  rQ  is  the  size  of  the  reflector  dish.  Equation  (27)  can  be  written 


F(s,^)exptjkQLa  +  jk^rsin©'  ((j)-*')  ]  sdsd<|>  (28) 


■I  (C  cos  n<b+D  sin  n*$)Pn(l-2s2)  (29) 

nm  nm  m 

n,m 

where  Pn(l-2s^)  is  a  modified  Jacobi  polynomial  which  relates  to  the  con- 
ventional  Jacobi  polynomial  by  the  following  relationship 

Pn(l-2s2)  -  /  2(n+2m+l)  P  (n,o)  (l-2s2)sn.  (30) 

m  m 


1  ^2ir 


I  =  kr 


U 


Let 


F(s,<J>)exp[jkoLa] 
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The  modified  Jacobi  polynomials  have  the  orthogonal  property,  that 

,1 


/: 


Pm  (1'2s)  *£.U-2s2)sds  =  5 


mm 


(31) 


Where  6  ,  is  the  delta  function.  The  Jacobi  polynomials  form  a  com¬ 

plete  set  similar  to  that  of  a  Fourier  series  which  also  forms  a  complete 
set.  These  functions  can  be  used  to  represent  any  arbitrary  function  if 
proper  coefficients  are  chosen.  Insert  this  relation  into  equation  (28), 
one  has 


2  c1  r2v 

I  =  kr  |  I  (  I  C 

o  I  I  n,m  n, 

Jo  Jq 


(cos  sin  n<$) 

m  nm 


(32) 


*  exp [  jk  r  s  sinfl 'cos  (<£-$')] d<j>  *  pn(l-2s2)  sds 
o  o  m 


By  use  of  the  identity, 

2tt 

cosnye 


L 


-jCcosp  2 

dp  =  —  j  (Z-) 

.n  n 
3 


(33) 


Equation  (32)  becomes 
1 

I  -  kr2 


2tt 

I  —  (C  cos  n<HD  sin  n$)  J  (k  r  s  sin©1) 
n,m  ^n  nm  nm  n  o  o 


.pn(l-2s2)sds 

m 


(34) 


By  use  of  the  identity 


■J1 


n+1/2  _  (n,o) 


y  1/2 


m 


n+2m+l 


(ay) 


2  1/2 
(l-2s  )  J  (yas) (ys)  ds 
n 


n  >  -(m+1) 


(35) 


Equation  (34)  becomes 

I  *  kr2  2*rr  I  I  ““  ^  2  (n+2m+l)  (C  cos  n$+D  sin  nt) 

o  . n  mn  m 

n  m  j 


Jn+2m+l(VoSine,) 


(36) 


k  r  sin0' 
o  o 
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The  C  and  D  are  evaluated  as  follows: 
mn  mn 


In  this  formulation,  the  integrand  I  is  evaluated  at  a  certain  0  angle. 
For  field  points  at  the  vicinity  of  this  angle,  this  integral  gives  a  very 
good  approximation.  However,  for  points  which  are  remote  from  this  angle, 
the  phase  error  becomes  intolerable.  To  reduce  this  error,  the  expotential 
term  of  Equation  (37)  may  be  expanded  into  a  Taylor  series.  When  C  and 
D  are  evaluated,  the  integration  can  be  carried  out  term  by  term.^^igh 
or§er  terms  can  be  calculated  by  use  of  an  iterative  process.  Hence,  the 
integration  needs  to  be  performed  only  once.  This  procedure  will  be  shown  as 
follows . 

For  a  parabloid. 


-  W«2  s2 

pcose  =  I  JT  “ 


Hence 


a  ro  (  2  4f“ 

exptjk^  =  exp^  jkQ (1-cos  6  )  If  \  ~  7? 

-  l  -  jH(s2-s2)  +  -4*r-  (s2-sV  • 


where 


H  “  ko  4?  (1_cos0  } 

e  2f 

s  ■  — 
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The  coefficient  C  and  D  takes  the  following  form 
mn  mn 


=  |  G(s,m,n)p”  (l-2s2)sds-j 


mn 


°  2  r1 

G(s,m,n)p”  (1- 


r1 

-jH  G(s,m,n)Pn 

J 


\l-2s  ) sds 


(l-2s  ) sds 


(40) 


where 


G  (s 


r1 

,m,n)  =  F(s, 
**  o 

_  r  i 


♦> 


cos  n<j> 
sin  n<j> 


d<|> 


Let  (B  )°  =  G(s,m#n)pn  (l-2s2)sds 
mn  m 


Then  (B  ) 
mn 


o 

x  r* 

*  1  G(s  ,m,n 


) (s2-s2)pn  (l-2s2)sds 
a  m 


(41) 


(42) 


From  Magnus  [7] ,  one  has 

.  2  2X  n  n  .  2.  n 

(s  -s  )P  =  a  P  +  (b  -s  )P 
am  mn  m-1  mn  a  m 


+  C  P  , 
mn  m+1 


(43) 


where 


m(n+m) _ 

mn  (n+2m) (n+2m+l) 


(44) 


b 

mn 


(m+n) 2 _ 

(2m+n) (n+2m+l) 


(44b) 


r  *  (m+1)  (n+m+1) 

inn  ~  (n+2m+l)  (n+2m+2) 


(44c) 
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Therefore 


and 


(B  )1  -  a  (B  )°  +  (bm  -  s  2  (B  )° 
mn  inn  n-i,n  mn  a  inn 


+  C  (B«4-1  n> 
mn 


(B  )P  -  a  (B  )P_1  +  (b  -  s  2)  (B  )p_1 
mn  mn  m— l ,  n  mn  a  mn 


(45) 


+  Cmn^Bmn+l,n)p_1 


(46) 


It  is  evident  we  may  use  this  iteration  process  to  compute  all  higher 

order  B  's.  Hence,  this  correction  can  be  carried  out  to  any  desired  degree, 
mn 


APPROXIMATE  SOLUTION 


J-(k  r  sin  9') 

Note  in  Equation  (36),  that  the  leading  term  (m*0,n*0)  is  — : — — - 

k.  r  sin  u 
o  o 

the  pattern  of  an  undistorted  uniformly  excited  equivalent  circular  aperture. 
All  the  remaining  terms  contribute  zero  at  sin0'*O.  The  higher  order  terms 
are  thus  perturbation  of  a  well-shaped  leading  term  beam.  This  term  there¬ 
fore  can  be  used  to  approximate  the  radiation  pattern  at  the  vicinity  of  the 
main  beam.  Since  off-focal  point  feed  shifts  the  beam  point  direction,  we 
shall  find  this  angle  deviation. 

Let's  assume  that 

6'  -  6  +  A  0  (47a) 

*’  -  *  +  A  *  (47b) 

If  A  6  and  A  *  are  small,  equation  (24a)  and  (24b)  can  be  approximated  by 
the  following  relations: 


A  0  cos  0  cos  *  -  A  *  sin  9  sin  *  *  e  /2f 

x 


(48a) 
(48b) 

Solving  A  0  and  A  *  from  these  two  simultaneous  equations  one  finds 

A  0  • 


A  9  cos  0  sin  $  +  A  *  sin  9  sin*  -  /2f 

y 


cos*  t  lit  +  sin*  e  /2f 

-x  I 


A  * 


cose 

cos*  e  /2f  -  sin*  e  /2f 
_ x _ y 

sine 


(49a) 

(49b) 
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These  two  equations  can  be  used  to  locate  the  shifted  antenna  pattern. 

COMPUTER  PROGRAM 

A  computer  program  which  computes  the  radiation  pattern  of  a  reflector 
disk  using  the  approach  described  in  this  report  is  compiled  in  FORTRAN 
language  and  has  been  run  on  an  ASC  machine  at  NRL.  Five  input  data  cards  to 
this  program  are  required  to  specify  the  parameters  of  the  dish  to  be  com¬ 
puted.  They  are  listed  as  follows: 

1.  The  first  card  specifies  the  axis  rotation  angles  of  the  coordinate 
of  the  feed  to  the  coordinate  of  the  reflector. 

2.  The  second  card  specifies  the  sequence  of  these  coordinate  rotation 
described  in  1. 

3.  The  third  card  specifies  the  displacement  of  the  feed  from  the  re¬ 
flector  focal  point. 

4.  The  fourth  card  specifies  the  number  of  terms  of  the  orthogonal 
series  to  be  used. 

5.  The  fifth  card  defines  the  size  of  the  disk  and  the  distance  of  the 
focal  point  and  also  specifies  the  feed  polarization. 

Figure  2  shows  an  example  of  such  a  computed  pattern.  This  figure  shows 
respectively  x,  y  and  z  components  of  the  radiation  of  a  dish  of  a  size, 

R  ■  12. 5A  and  f  *  12. 5A.  The  feed  horn  is  linearly  polarized  at  Y  -  direction 
and  has  a  taper  of  10  dB.  Location  of  this  feed  is  displaced  from  the  focal 
point  in  x  direction  by  an  amount  of  £  -  3.18A.  KAU  *  2  x  ir  xR  x  sin0. 

Where  0  is  the  angle  of  a  field  point.  The  plot  is  at  a  cut  of  4>  *  0. 

To  achieve  this  plot  a  total  of  200  field  points  has  been  computed. 

Total  computer  time  of  this  plot  including  compiling  and  evaluating  coef¬ 
ficient  is  41  seconds  on  an  ASC  machine. 
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